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ABSTRACT 


The problem considered Is the design of a feedback 
system containing a linear, time-invariant, minimum-phase 
plant, whose parameters are known only within given bounds, 
such that the time-response of the system remains within 
specified limits. A quasi-opt imal design, for given design 
constraints, is one which minimizes the effect of white 
sensor-noise on the input to the plant. Horowitz and Sidi 
have presented a method for the design of such a quasi- 
optimal, linear time-invariant system. 

This report investigates the use of the non-linear device 
known as the Clegg Integrator in the design of such a system. 
The describing- function of the Clegg Integrator has the same 
magnitude characteristic, apart from a scale factor, as the 
linear integrator, but has 52° less phase-lag, at all fre- 
quencies, than the linear integrator; thus, when used in a 
feedback system, it provides a larger stability margin than 
the linear integrator. This property allows the non-linear 
feedback system to be designed so that the sensor-noise is 
attenuated more than in the linear design. 

The non-linear time-response is calculated in terms of an 
equivalent linear representation, and the reduced phase-lag in 
the frequency-domain is seen to correspond, in general, to 
smaller overshoot in the time-response. In order to obtain 
satisfactory steady-state response, the non-linear device is 
modified by putting a linear integrator in parallel with it. 
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It is shown that for given bounds on plant-ignorance and 
time-response tolerances, the non-linear feedback system can 
be designed such that the effect of white sensor-noise on the 
plant-input is smaller than in a linear design. A method is 
developed for implementing the non-linear design in terms of 
an equivalent linear design by Horowitz and Sidi’s method. 



CHAPTER I 


STATEMENT OP THE PROBLEM 


1.1 Introduction 

An Important, fundamental problem In control Is the design 
of a system whose performance remains within specified limits 
despite the uncertainty in the parameters of a vital system 
portion denoted as the plant. In Pig. 1.1, let P(s) be the 
transfer-function of a linear, time-invariant, minimum-phase 
plant whose parameters assume values within a given range. 

Note that we are not considering the case of a time-varying 
plant whose parameters change during the time an input signal 
is present; rather, we suppose the plant to be constant while 
a signal is being applied, but allow it to change between 
successive applications of the input; or else, we may take 
the plant parameters to be fixed but known only within certain 
tolerances. 

Given such a plant P(s) with the range of values of its 
parameters, and given a range of acceptable overall system 
performance, such as, for example, bounds for its step-response 
as in Pig. 1.2, the problem (for design with linear, time- 
invariant elements) is to find fixed compensating functions 
G(s) and P(s), such that the response of the closed-loop system 
of Fig. 1.1 lies within the specified bounds for all values 
of the plant-parameters within the given range. 
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1 . 2 Principle of Linear Design^ ^ 

The general principle underlying a linear design for 
this problem can be described as follows: the closed-loop 

transfer- function in Pig. 1.1 is 


A 

T(s) = 


C(s) 

R(s) 


N( s ) =0 


F(s) L(s ) 
1+L( s) 


where 


L(s) 


A 


P(s) 


G(s) 


With s = jo), if | L( j to) | >> 1, then T(jco) ^ FQq) , and is thus 
quite insensitive to changes in L(joi) (caused by changes in 
P(jto)), so long as |L(jw)| remains much larger than 1. In a 
linear design, a minimum-phase G(s) is so chosen that 
|L(j w )| >> 1 for all possible plant-conditions, over a 
frequency-range 0 < to < to whose width to depends on the 
response specifications. Thus T(jto) % F(joo) for 0 < to < to c , 
independent of the actual plant condition, and if to is 
properly chosen, the closed-loop time-response essentially 
depends on F(jco) and can be made quite insensitive to changes 
in the plant. 

We note from Fig. 1.1, however, that the transfer- 
function between the noise-source n(t) in the feedback sensor 
and the input z(t) to the plant is 


Z(s) 

1 

L(s) 

N(s) 

R( s ) =0 

_1+L(s)_ 


B ( s ) 
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If 


S nn M 

s zz <“> 

s zz <“> 


rms 


= spectral density of n(t) 

= spectral density of z(t), then, 
= s nn ( w ) I B( jw) | 2 , and 

1/2 


r- 

r oo _. 

1_ 

S (w)do) 
zz 

— GO __ 

2tt 

— j 


Considering B(jco) for m >> u) , where | L ( J oj) J << 1, 

|B(jw)| 2 % | 2 = |G(jo3)| 2 , and if |G(jto)| > 1, the sensor 

P(jo>) 

noise-components at such frequencies are amplified at the Input 

to the plant and tend to produce saturation. 

| G ( J u) ) | > 1 implies that | L ( j cc ) | > |P(jw.)|, a condition 

that occurs when the plant is unable to provide the required 

loop-gain which must then be supplied by the compensation 

G(jm). The condition | L ( J u) ) j > |P(jw)| can occur over a wide 

frequency-range even when |L(ju>)| < 1 over that range, as 

illustrated in Fig. 1.3 with P(s) = — and 

s 


L(s) = P(s)G(s) ■ 2048 


(s+1) 


s(s+8) (sfl6) 2 


| L ( J oj) j < 1 for all w > 8 rad/sec but |G(j<o)| > 1 until co 
exceeds about 40 rad/sec. Fig. 1.4 shows |G(jw)| as a 
function of u on linear scales, and it is seen that the 
frequency-range 8 < w < 40 (throughout which L(jto) < 1) is, 
in fact, the significant range for the amplification of sensor- 
noise at the plant-input. 

It is, therefore, desirable to reduce | L ( j to ) | as rapidly 
as possible beyond the frequency-range 0 < w < w , in which 
sensitivity-specifications require it to be; large. However, 
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in a minimum-phase transfer-function, rapid attenuation in 
the magnitude characteristic is accompanied by a large phase- 
lag, and a limit is set to the rate of attenuation of | L ( J to) ( 
by the stability requirement that the phase-lag of L(jw) cannot 
exceed 180° till J L ( J to ) | has become less than 1, for all 
possible plant-parameter values. Thus, as a result of the 
relation between the magnitude and phase functions of L(jw), 
the attenuation of ( L ( j to ) j has to be gradual until j L ( J to ) ( has 
become less than 1 at a phase-lag less than 180° , for all 
plant-conditions. 

We may define an optimal design L(s) for the above 
problem as one which satisfies the specifications on closed- 
loop response while minimizing the effect of sensor noise on 
the input to the plant. This definition takes into account 
the power spectrum of the noise source. This problem has not 
yet been solved. Another problem, which approximates the 
above when the sensor-noise is white, is defined as follows: 

Let q = (number of finite poles - number of finite zeros) of L(s). 

Then, for large |s|, L(s) has the asymptotic form 

L(s) % (1.1) 

s q 

For a given q, the optimal L(s) is defined as that which meets 

the specifications on closed-loop response and has the smallest 

possible value of k in Equation (1.1). 

This problem of optimal linear design has been investi- 

T 1 2 "1 

gated by Horowitz and Sidi ’ , who have established the 

existence and properties of the optimal design for a class of 
problems, and developed a procedure for its implementation. 



1. 3 Design with Non-Linear Elements 

The fixed relation that exists between the magnitude and 
phase of a transfer- function sets a limit to the maximum rate 
of attenuation of | L( Jui ) | possible in a stable linear 
system, over a frequency-range that is important in regard to 
amplification of sensor-noise. Can this limitation, inherent 
in a linear design, be overcome by the use of non-linear 
devices in the compensation? Suppose there is a non-linear 
device with a describing-function in which the magnitude 
characteristic is associated with a smaller phase-lag than 
that due to the same magnitude characteristic in a linear, 
time-invariant device. Such a device would appear to permit 
a sharper reduction in the magnitude of the loop describing 
function than that allowed in a linear design. We are then 
led to the question, "Using such a non-linear device, can the 
specifications on response, in the problem stated above, be 
met with a compensation which produces a smaller noise output 
than the optimal linear design for the same problem?" 

This report presents the results of an investigation of 

this question, in which the non-linear device is the Clegg 
T3 41 

Integrator It is shown that the non-linear device 

leads to a design with a smaller effect of sensor noise than 
in the optimal linear design, for problems of two-degree-of- 
freedom design. 




CHAPTER II 


ANALYSIS OF NON-LINEAR RESPONSE 


2 . 1 The Clegg Integrator 

n 41 

This device 1 - 5 is a non-linear integrator obtained 
from a linear integrator merely by arranging its output to 
be reset to zero whenever the input signal crosses zero and 
changes sign. Between consecutive zero-crossings of the 
input, the device is just a linear integrator. It can be 
simulated on the analogue computer by using logic circuits 
to detect zero-crossings of the input and to switch a linear 
integrator between "reset" and "compute" modes. The simula- 
tion we used is shown in Pig. 2.1. 

Describing Function of the Clegg Integrator 

If a non-linear device produces a periodic output y(t) 
to the sinusoidal input x(t) = A sinwt, its describing- 
function is defined as the complex ratio 



A 

where B sin(wt+0) is the fundamental component in the Fourier 
series expansion of y(t). 

When the input to the Clegg Integrator is x(t) = A sinwt, 
its output y(t) has the periodic waveform shown in Fig. 2.2, 
and is given, over the first period, by the expression 


8 
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< A . , A 

- — COSCOt + — 


0) 


CO 


y(t) =< 


A A 

- - C OSCOt - - 
CO (0 


0 < cot < TT 


TT . < COt < 27T 


i.e. y(t) = . - ^ coscot + f(t) , where f(t) is a square wave 
A 2ir 

of amplitude — and period — . 

<o to 

The Fourier series for f(t) is 
4 a 


f(t) = 


TT CO 


sincot + sin3<ot + — sin5cot + . . 
3 5 


and, therefore, y(t) has the Fourier series 


y(t) = 



r, — i 

4 a a 

— sincot - — coscot 

+ M 

- sin3wt + i sin5wt + ... 


irco co 

TTCO 

L3 5 J 


The fundamental component of y(t) is 

1/2 


[4 , , 1 

A 

, 16 

— sincot - coscot 


1 + —■ 

TT 

co 

TT 2 


sin cot - tan 


-1 


= 1-62A -j 38° 


co 

S Be’’ 9 


Therefore, the describing function D(jco) of the Clegg Integrator 
is given by 


D(jco) 


b e j e _ 1^62 e -j 38 ° 

A w 
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while the transfer function of a linear integrator, for 

s = jo), is ( — ) . We note that the describing function of 
\ ju>/ 

the non-linear integrator has the same magnitude character- 
istic as the linear integrator, apart from a scale factor of 
1 . 62 , but has a lag of 38° at all frequencies, which is a 
lead of 52° over the phase of the linear integrator. The 
describing-function is independent of the input amplitude A, 
a fact which makes closed-loop calculations as simple as 
open-loop calculations. 


2.2 Stability in the Absence of Input 

Prom the describing-function of the Clegg Integrator, 

it appears that a closed-loop system that Is unstable with a 

linear integrator might be stable when a Clegg Integrator 

replaces the linear integrator. For example. In the linear 

1c 

system shown In Fig. 2.3-a, the loop transmission L(s) = — 

s 

and the condition for free oscillations Is 1 + L(jw) = 0, 

1/2 

whose solution is w = k ; hence, the system is unstable 
for all k, the frequency of oscillations depending on k. 

The result of analogue simulation for k = 1 is shown in 
Fig. 2 . 3-b . 

The system of Fig. 2.4-a is obtained by replacing a 
linear Integrator In the system of Fig. 2.3-a by the Clegg 
Integrator. Since the loop-describing function 




1. 62ke 


J52° 


( jw) 


, the equation 1 + L ri (ja>) = 0 Is not 


satisfied by any real w and k. Thus, according to 



describing- function theory, in the absence of external input 

the non-linear system should be stable for all k. This is 

borne out by the results of analogue simulation for k = 10 

and k = 100, shown in Figures 2 . 4-b and 2.4-c, respectively. 

As another example, consider the linear system of 

Fig. 2.5-a. The loop-transmission L(s) = — — - and the 

s(s+l) 2 

equation 1 + 'L(jw) = 0 has the solution to = 1, k = 2. Thus, 
the system Is stable for 0 < k < 2 and breaks into unstable 
oscillations of frequency 1 rad/sec at k = 2. The result 
of analogue simulation for k = 2.1 is shown in Fig. 2.5-b. 

On replacing the linear integrator in the system of 
Fig. 2.5-a by the Clegg Integrator, the non-linear system 
of Fig. 2.6-a Is obtained. The loop-describing function 

1 62ke J52 ° 

L,(joo) = — — - , and the equation 1 + L.(jto) « 0 

d j<o(jw+l) 2 d 

has the solution to = 2.9 and k = 17. The non-linear system 

is thus stable for 0 < k < 17 and breaks into unstable 

oscillations of frequency 2.9 rad/sec at k = 17. Figures 

2 . 6-b and 2.6-c show the non-linear system responses at 

k = 5 and k = 1 6 respectively, and there is evidence of 

a stable limit cycle, with an amplitude of about 0.15 volts 

for the case k = 16, while Fig. 2.7 shows the unstable 

oscillations at k = 18. 

The stable limit cycle occurring in Fig. 2.6-c, for 
example, is not predicted by the single input describing- 
function of the Clegg Integrator. It occurs in a situation 
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where the oscillations at the input to the Clegg Integrator 
have a non-zero mean value and a dual-input describing- 
function of the Clegg Integrator must be investigated to 
account for it. This point is considered again in Sections 
2.3 and 2.4. 

2 . 3 Equations for Calculation of Forced Response 

The verification of stability in the absence of input, 
however, gives no indication of the forced-response of the 
non-linear system, such as, for example, its response to a 
step-input. We now consider the calculation of forced- 
response of a system containing the Clegg Integrator, assuming 
that all other elements of the system are linear and time- 
invariant . 

We have the following relation between the input x(t) 
to the Clegg Integrator and its output y(t) : 


y(t) 


r t 


x(p)dp 


( 2 . 1 ) 


where t_ is the instant of the last zero-crossing of the input 
z _ 

x prior to t. It is thus clear that the zero-crossings of the 
input have a decided effect on the output. Consider now the 
interval [t Q , t^J and suppose that t^,..., t n are the instants 
of zero-crossing of x(t) in this interval, with t Q < t-^ < tg-.. < 
t n < t f . Then, over the sub-interval t^ < t < t^ + ^. 
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(i = 0,1,.. . ,n-l) , we have 


y ( t) = 


X ( p ) dp 


- 42 “ * 1+1 


t . 

X 


r t- 


x(p)dp - 


x(p)dp 


) t 


o 


x ( p ) dp - 


o 


V 

k=l J t 


x ( p ) dp 


k-1 


and, = 0, the output being reset to zero as x(t) 

crosses zero at Using the unit-step function u(t), we 

may write a single expression for y(t) for the entire interval 
[t Q , t^] as follows: 


^ t 


y (t) = 


n 


x(p)dp - 


k=l 


x(p)dp 


'k-1 


u(t-t k ); t o <t<t f (2.2) 


We note that y(t) is discontinuous at the instants t i of zero- 
crossing of x(t) , being reset to zero at t^ from the value 

t . 

1 A 

x(p) dp = (2.3) 

' t i-l 


y( t i -). 


Substituting from (2.3) into (2.2), 
t n 


y(t) = 


x 


( P ) dp Y u(t-t 1 ) 


i=l 

- 1 i — n 


x(p) - Y 


i=l 


dp , t „ < t < t . 

o - — ] 


(2.4) 





which can be represented as the output of a linear integrator 

to which is applied the sum of x(t) and an impulse-train 

a(t) = - £ y^Mt-t.). We thus arrive at the equivalence shown 
i 


in Pig. 2.8, in which the Clegg Integrator is replaced by a 
linear integrator and the resets are represented by the addi- 
tion of a suitable impulse-train at its input. 

With this linear representation of the Clegg Integrator, 
we may calculate the response as the sum of the linear responses 
to the actual input and to the "equivalent input" consisting, 
of a train of impulses. 

Thus the closed-loop non-linear system of Pig. 2.9 may 
be represented by the linear system of Fig. 2.10, and the 
forced-response calculated by adding the response due to the 
input r(t) to that due to the impulse-train. We thus obtain 
the pair of equations 


c(t) 
y( t) 


^ t 

> o 

" t 

> o 


r (p)h( t-p) dp - 


r(p)v(t-p)dp - 



(2.5) 

( 2 . 6 ) 


where 


H ( s ) = 

fh(t) - C(s) 

_ C( s ) 

« 

3- R(s) 

A(s) 

V(s) = 

Cv(t) = Y(s) 

_ Y(s) 

<0 

L- R(s) 

A(s) 


for the linear 

> 

system of Pig. 2.10; 

/ 


and. 



lim 

t+t^- 


y(t) , 


with the t^'s being the instants at 
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which e(t) = 0, i.e. r(t^) = c(t^), i = 1,2,,..; these are 
the instants at which the Clegg Integrator is reset. We 
note that for t < t^, hCt-t^) = vCt-t^) = 0. 

Equations (2.5) and (2.6) do not seem to admit of a 
general analytical solution for c(t) and y(t) but a numerical 
solution may always be constructed in the following manner: 
for 0 < t < t x (i.e. till the first reset), c(t) and y(t) are 
given by just the integral terms in Equations (2.5) and (2.6): 


c(t) 
y (t) 


f t 

= r (p )h( t-p ) dp 
> o 


t 

= r (p ) v( t-p ) dp 

J o 


(2.7) ) 


> 0 < t < t^ 

( 2 . 8 ) 


Expression (2.7) for c(t), in fact, determines whether, at 
all, there is an instant t^> 0 at which ett.^) = r(t 1 )-c(t 1 ) = 0 
If that condition never occurs, then the input to the Clegg 
Integrator never changes sign and it Is never reset; in that 
case, c(t) and y(t) are identical with the corresponding linear 
responses (those obtained when a linear integrator is used in 
place of the Clegg Integrator) , and are given by Equations 

(2.7) and (2.8) for all t > 0. If, however. Equation (2.7) 
does give a t^ > 0 at which r(t^) - c(t^) = 0, then Equation 

(2.8) determines the value of y^ = y(t^-), and the terms 
-y 1 h(t-t 1 ) and -y 1 v(t-t^) are added to the right-hand side of 
Equations (2.7) and (2.8) respectively, enabling the non-linear 
response to be calculated for t > t^. 

In general. Equation (2.5) determines whether there is a 
reset at t^ +1 following a reset at ty If no reset occurs 
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If another reset does occur at tj + ^, then Equation (2.6) 
determines -yj + 1 = -y(tj +1 -) , the strength of the impulse to 
be added at the input to the integrator at tj +1 in order to 
calculate c(t) and y(t) for t > tj + ^. 

Setting r(t) = u(t) in Equations (2.5) and (2.6), we get 
equations for the step-response of the system of Pig. 2.9. It 
is obvious that if the step-response of the corresponding 
linear system (obtained upon replacing the Clegg Integrator 
by a linear integrator) has no overshoot, then the condition 
c(t) = r(t) = 1 never occurs for finite t; thus t 1 does not 
exist and the non-linear step-response is identical with the 
corresponding linear step-response. 

What is the nature of the non-linear step-response in the 
case when the corresponding linear step-response has overshoot? 
The instant of the first overshoot in the linear step-response 
is the instant t-^ of the first reset of the Clegg Integrator 
in the non-linear case, and for 0 < t < t^, the linear and 
non-linear responses are identical. Equations (2.5) and (2.6) 
enable the further calculation of the non-linear response. 

It is instructive to carry out the calculations for a second- 
order system in which the corresponding linear step-response 
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can have overshoot. An analytical solution of Equations 

(2.5) and (2.6) is possible in this special case. 

Let us consider the step-response of the non-linear 

system shown in Fig. 2.11. The corresponding linear system, 

C(s') 

shown in Fig. 2.12, has the transfer- function H(s) = = 

— P , and will have overshoot in its step-response if 

s^+as+k 

a 2 

k > — . The linear step-response is given by 


c(t) = 1 - — e ^n 1 ^ S in 


/ 1-6 2 


(^ n / 1-6 2 t + cos" 1 ^ (2.9) 


where 


"n = / ~ i ‘ 


6 = 


a 


2/~k 


a 
2 uo 


n 


the instant of first overshoot, which is the instant t-^ of 
the first reset of the Clegg Integrator in the non-linear case 
is given by 


h = 


7r-cos - ^<5 


(2.10) 


to 


n 


Zl-6 2 


On calculating the expressions in Equations (2.5) and (2.6), 
it is found that 


t , - 1, + 

1 /T& 


to - 

n 

iTt-cos - ^ 


to 


n 




, i = 1,2,3,... 


, i — 1,2,3,..* 


( 2 . 11 ) 
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and 

c(t) = 


n 

1 - sin 


/l-6 2 




t + cos~^S ) » 0 < t < t- 


1 - sin 

/l-<5 2 


V> n /l-& 2 (t-t 1 ) 


, t, < t < t, 


C t - 


co 




n 


,/l-6 2 


3 t > t, 


( 2 . 12 ) 


Equation (2.12) describes a step-response with a periodic 
undershoot, and the result of analogue simulation shown in 
Fig. 2.13 confirms this conclusion. 

We thus find that while the non-linear system of Fig. 2.11 
is stable without an input, as one expects from the describing- 
function of the loop, its response to a step-input contains 
forced oscillations and does not reach a steady state. This 
result is not explained by the single-input sinusoidal 
describing function of the Clegg Integrator for the simple 
reason that it is inadequate for the situation in which a step- 
input is applied. The actual input to the Clegg Integrator is 
e(t) = 1 - c(t) and it is seen to oscillate about a non-zero 
mean value. The sinusoidal describing-function is not 
applicable in this situation. 

We have shown earlier that the non-linear step-response 
is identical with the corresponding linear step-response 
whenever the latter has no overshoot. In the second-order 
example above, we also find that if the linear step-response 
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does have overshoot, the non-linear response contains sustained 
oscillations. Simulation of several higher-order systems 
suggests that, in general, the non-linear step-response con- 
tains sustained oscillations whenever the corresponding linear 
step-response has overshoot. An analytical proof of this fact 
has not been found, but the instances in which it has been 
verified are themselves sufficient to show that the Clegg 
Integrator cannot be used without some modification. 

2 . 4 A Modified Non-Linear Device 

It has been found that when a linear integrator is put in 
parallel with the Clegg Integrator, the sustained oscillations 
in the step-response are eliminated. The modified non-linear 
device, shown in Pig. 2.l4-a, has the describing-function 

^ , which has less phase-lag than the 

corresponding linear transfer- function f — 

\ 

The phase lead obtained In the non-linear device over the' 
corresponding linear transfer-function is plotted as a function 
of b in Fig. 2.l4-b. 

Step-Response 

Fig. 2.15 shows a closed-loop system with this non-linear 
device and Fig, 2.16 shows its step-response. The corresponding 
linear system Is shown in Fig. 2.17 and its step-response in 
Fig. 2.18. The non-linear response now attains the same steady- 
state as the linear response but with much less oscillation 
than the linear response. By simulation of a third-order 



( 1 . 62 e J52 °+b 
\ Jw 
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system, it was also verified that the limit of stable operation 
of the non-linear system corresponds to the sinusoidal 
describing- function ^ — ~ j of the non-linear device, 

both in the absence of input and the presence of a step-input . 

In explanation of this result, we note that due to the linear 
integrator in parallel with the Clegg Integrator, any sustained 
oscillations at the input to the Clegg Integrator must have 
zero mean, whether or not there is an input (unless the response 
itself is unbounded) . The sinusoidal describing- function is 
applicable in this situation and the results agree with con- 
clusions based on it. In other words, within the stability 
limit set by the describing- function, the step-response cannot 
have sustained oscillations. 

Ramp-Response 

While the presence of a linear integrator in parallel 
with the Clegg Integrator inside a closed-loop removes sus- 
tained oscillations from the step-response of the closed-loop, 
the nature of the ramp-response remains to be considered. Of 

course, with the device (Clegg Integrator + — ) at the error- 

s 

junction, if the corresponding linear ramp-response has no 
overshoot, then the non-linear ramp-response Is identical with 
the linear ramp-response and is thus free of sustained oscil- 
lations. However, if the corresponding linear loop-transmission 
is of Type-0 or Type-1 and the corresponding linear ramp- 
response has overshoot, it turns out that the non-linear 
ramp-response contains sustained oscillations, as demonstrated 
by Pig. 2.20, the ramp-response of the system of Pig. 2.19. 
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This happens because , with a ramp Input, oscillations at the 
Input to the Clegg Integrator do not have zero mean-value If 
the loop-transmission is of Type-0 or Type-1, but must have 
a non-zero mean value to give rise to the ramp in the output. 
The sinusoidal describing-function of the Clegg Integrator is 
inapplicable in this situation. Following this train of 
thought, we are led to conclude that if the loop-transmission 
is of Type-2 (or higher) the ramp-response should be free of 
sustained oscillations, and, in general, if a steady-state 
condition requires that at large t, the input to the Clegg 
Integrator has zero mean, then the sinusoidal describing- 
function is applicable, and the non-linear system response 
will be free of sustained oscillations within the stability 
limit determined by describing- function theory. This conclu- 
sion is borne out by Fig. 2.22, the ramp-response of Fig. 2.21. 
We note that since the linear system corresponding to the 
system of Fig. 2.21 is unstable, the corresponding linear 
ramp-response is, in fact, unstable, while the non-linear 
response is stable. 

2 . 5 Equations for Calculation of Forced-Response with 

Modified Non-Linear Device 

The calculation of time-response with the modified non- 
linear device follows directly from the linear representation 
suggested earlier for the Clegg Integrator, and is illustrated 
with the system of Fig. 2.23. The equivalent linear repre- 
sentation is given in Fig. 2.2k. The transfer-functions needed 
in the calculations are: 



(i) 


(ii) 


(iii) 


(iv) 


H(s) 


V(s) 


H a (s) 


V s) 


£ h( t) ft 
£v(t> ft 

=£h a <t) 
= £ a (t> 


CUi = _L(s)_ } wlth L(s) = P(_g.)ll±b) 

R(s) 1+L(s) s 

Y(s) = 1 f 1 \ 

R(s) s \l+L(s ) ) 

4 °( s ) = 1 l L(s) \ = H(s) _ C h(t ) 

A(s) l+b\l+L(s )/ 1+b J 1+b 

= VCs) + b 

1+b s ( 1+b ) 


A Y ( s ) = 1 
A(s) s 


1 + 


P(s)b 


1+L(s) 


v( t)+b 
1+b 


The following equations, derived In the manner of Equations 
(2.5) and (2.6), describe the time-response In this case: 

r t 


c(t) 


r(p)h(t-p)dp 


y y 1 K a (t-t 1 ) 


(2.13a) 


J o 


1 


y (t) 


( ^ 

r(p)h(t-p)dp - 

Y y i h(t - t i ) 

(2.13b) 

o 


1+b ^ 1 1 


r t 

r(p)v( t-p)dp - 

Y y i v a (t “ t l ) 

( 2 . 14a) 

0 


1 


r t 

r(p)v(t-p)dp - 

Y y, [b+v( t-t . )1 

( 2 . l4b ) 

0 


1+b ~ 1 L 



where, as before, 

y. = 11m y ( t ) 
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Equations (2.13a) and (2.14a) differ from Equations (2.5) 
and (2.6) only in the appearance of h„(t) and v (t), in place 

a. a. 

of h(t) and v(t), in the summation on the right-hand side. 
When b = 0, of course, Equations (2.13) and (2.14) reduce to 
Equations (2.5) and (2.6), since, in that case, v(t) = v (t) 
and h( t) = h ( t) . 

3 


2 . 6 Second-Order Non-Linear Step- and Disturbance-Response 

Applying Equations (2.13) and (2.14) to the second-order 
case shown in Pig. 2.15, we now derive an expression for the 
maximum overshoot in the step-response, which occurs between 
t-^ and t 2 , the first and second resets of the Clegg Integrator. 


With k e = k( 1+b ) , 
H(s) ‘A 

R(s) 


k 

e 

s^+as+k 

e 


and 


h( t) 


/ 1-6 2 


e - w n 5t 


sin 



where w = /k~ 
n e 


6 = 


2/k~ 

e 


-lx 

TT-COS O 

to / 1-5 2 
n 


(2.15) 



Using Equations (2 . 13) , (2.14) and (2.15), the step-response, 
for t^ < t < t 2 , is given by 

6 


c ( t ) = 1 - 


b-26e 


/ 1-6 2 


(ir-cos - ^6 ) 


( 1+b ) /l-5 1 2 


—to 6 1 . i 

e n sin to 


/l-6 2 t+cos -1 6 


n 


t, < t < tj 


cl c 

Setting — = 0 , we find that the maximum of c(t) occurs at 
dt 


max 


A -& 2 


, and 


to 


n 


0(t max> * 1 + 


b-2Se 


/l-6 2 


( Tr-cos _1 6 ) 


1+b 


-rr 6 


/l-S 2 


Thus, M n (6,b) = the maximum overshoot in the non-linear 

step-response 


— — (Tr-cos*"^ )~ ' 


b-2Ae 


/l-S 2 


1 + b 


7T 6 


JIT 2 


(2.17) 


Since M l ( 6) = the overshoot in the corresponding linear 
step-response 

tt<5 


= e 


/l-<$ 2 


(2.18) 


we note that for b > 0, M (6,b) < M T (6). Expression (2.17) 

n ij 

can be negative in certain cases, showing that in those cases, 
the non-linear step-response has undershoot instead of overshoot. 
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3M (6,b) 

In fact, < 0 for 0 < 6 < 1 (see Appendix) and, 

36 

therefore, for fixed b, the maximum undershoot occurs when 
6=1 (critical damping), and is given by 


Maximum undershoot = - M n (l,b) 


2e~ X = 0.736 
1+b l+b 


(2.19) 


For later use in non-linear design, we collect the 
following results on the second-order, non-linear step-response 


of the one-degree-of- freedom structure of, Fig. 2.15 


(i) Maximum overshoot = 


b-26e 


(tt-cos -X 6) 

/l-6 2 


l + b 


tt6 


fa 


D 7 *3 f\ 

(ii) Maximum undershoot (occurs when 6 = 1) = ■ 

l+b 


(2.17) 


(2.19) 


The 6 here, of course, refers to the damping in the corresponding 
linear transfer-function, and the maximum undershoot in the 
non-linear response occurs when the corresponding linear response 
is critically damped. 

It may also be noted that expressions (2.17) and (2.19) 
also apply to the disturbance-response of the one-degree-of- 
freedom structure of Fig. 2.15, since, for this structure 

disturbance-response = 1 - (step-response) . 

Equations (2.17) and (2.18) compare the overshoots in 
the disturbance-responses of the corresponding non-linear and 
linear systems of Figures 2,15 and 2.17. We note that for a 



given damping in the linear system, the non-linear response 
has less overshoot than the corresponding linear response; in 
consequence, for a given overshoot in the disturbance-response , 
the non-linear system can be designed with a smaller damping 
(in the corresponding linear system) than is possible in a 
linear design . 

This is the principle that is applied later to show that 
in a two-degree-of-freedom design to meet specifications on 
step-response and disturbance-response, a non-linear design 
is possible which has a smaller bandwidth of (the corresponding 
linear) loop-transmission and transmits less feedback sensor- 
noise to the plant than the optimal linear design of Horowitz 
and Sidi. 

In order to establish this result, we next consider the 
noise-response of the Clegg Integrator. 


2.7 Noise-Response of the Clegg Integrator and the Modified 
Non-Linear Device 

An input x(t) to the Clegg Integrator produces the output 


y ( t ) given by 


y(t) = 


x ( p ) dp 


( 2 . 1 ) 


where t is the instant of the last zero-crossing of x prior 
z 

to t. Let T(.t) = t - t 

Ci 

= the interval between t and the last zero- 
crossing of x prior to t; 

hence, t = t - T(t), and thus 

z 
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^ t 

y(t) =. x(p)dp 

J t-T(t) 


When x(t)- is a stationary random process, the interval 
T(t) is itself a stationary random process (or, for fixed t, 
a random variable). In principle, the statistics of T may 
be determined from those of x(t), but such calculations have 

rsi 

been possible only for a few processes , owing to the 
mathematical difficulties that arise. Even if the statistics 
of T are known, calculating the statistics of the output y(t) 
is a formidable task, since this requires finding the expected 
value of a stochastic integral, one of whose limits is a 
random process whose statistics, in turn, depend on the inte- 
grated process. 

F51 

It is known , however, that for a zero-mean Gaussian 

Markov process x(t), the mean and variance of T tend to zero 

as the process x(t) becomes more and more uncorrelated. This 

suggests an approximation that simplifies the problem and 

enables us to calculate the power spectrum of the output. In the 

ft 

expression x(p)dp, the random variable T in the lower 

t-T 


limit is replaced by its mean value T . We then have the 
approximate relation 


r t 


y(t) = 


x ( p ) dp 


( 2 . 20 ) 


t-T 


m 


which is a stochastic integral with deterministic limits. 



(i) Let x and y be the expected values of the 
m m c 

stationary processes x(t) and y(t) respectively. Using 
to denote expected value, we have, from Equation (2.20) 


y m = E(y(t) } 


= E { 


x ( p ) dp } 


t-T 


m 


r t 


t-T 


E{x(p) }dp 


m 


= T x 
m m 


in particular, y^ = 0 if x m = 0. 


(ii) Let F denote Fourier Transform and let 
R (x) = auto-correlation of x(t) 

XX 

R (t) = cross-correlation of y(t) and x(t) 
yx 

R (x) = auto-correlation of y(t) 

«y «/ 

S (w) = spectral density of x(t) = F(R yv (t)} 
xx xx 

S (w) = cross-spectral density of y(t) and x(t) = 

yx 

S (u>) = spectral density of y(t) = F{R (t)} 

«y *y o *j 


We have 


R (t) = E{y(t+x)x(t) } 
y A 

r t+x 


= E{ 


x(p)x(t)dp> 
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E 


F{R yx (x)} 


t+x-T 

m 
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r t+T 


t+x-T 


E(x(p)x( t ) }dp 


m 


r t + x 


R xx (p - t)dp 


t + x-T 


m 


m 


R xx (T - T m +v)dv 


(2. 


Taking Fourier Transforms on both sides, 


r 


S (u>) = 
yx v 


m 


J o 
r T m 


P(R xx (T -V v)1<lv 


3 U (,)^<-.) 4V 


* S xx M 


T 


m 


( jco(v-T m ) dv 


S (w) 

XX 


l__ e -jwT m 1 




( 2 , 


Then, using the technique of statistical linearization^ 5 ^ , 
we define G(jw), the random-describing function of the Clegg 
Integrator, as follows: 


A S v -x (u) 
G(ju) = -22- 


S xx<“> 




( 2 , 


r C n "1 

It then follows 1 5 J that, for the Clegg Integrator, 


21a) 


21b) 


21c) 
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S (to) ? m p / coT 

= 1 G( Jto) | 2 = i- sin 2 ( — 

S xx (to) to 2 V 2 


( 2 . 2.1d) 


For a linear integrator, the corresponding relation is 


s .vv (M) ■ L 
S xx<“> 


Using the relation 


rms 


— * 

f 00 — 

1_ 

S „ (oj) doi 

- 2Tr . 

yy 

_oo ___ 


1/2 


( 2 . 22 ) 


Equations (2.21d) and (2.22) enable us to conclude that the 

Clegg Integrator in cascade with a pure attenuation of 0.5 

gives a smaller noise output than a linear integrator — , 

s 

regardless of the actual value of T m , so long as Equation 
(2.20) constitutes a reasonable approximation to Equation (2.1). 

Turning to the modified non-linear device of Fig. 2.14, 
the input x(t) and the output y(t) are related by 


y ( t ) = b 


x(p)dp + 


x(p)dp 


t-T 


(2.23) 


where, as before, (t-T) is the instant of the last zero- 
crossing of x before t. Approximating T by T , as before. 


r t 


y(t) = b 


x(p)dp + 


x(p)dp 


(2.24) 


t-T. 


m 


r 


Writing y(t) = y 1 (t> + y 2 (t), where y 1 (t) = b 


x (p ) dp = 


= output of the linear integrator — , 

s 
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/ t 


y 2 ( t ) = 


x(p)dp = output of the Clegg Integrator, 


t-T 
. m 


we obtain. 


s yx<“) = s y lX ( “> + s y 2 x (w) 


= S ( to ) 
xx 


— + G(ju>) 


G(jto) = 


l- e ~ jtoT ni 

jw 


, where 


A S (to) 
Thus, G, ( j to ) = 

J- o 


S xx<“> 


= describing-function of the modified non-linear device 


b + l-e 


-JwT m 


jw 


( 2 . 25 a) 


and 


S yy^ 

s xx (w) 


= iG.Uw) 


h d + 4(l+b)sin £ 


u)T r 


( 2 . 25 b) 


CO 


For the corresponding linear transfer-function , we have 


S yy^ _ (b + 1) 2 _ b 2 +( 1+b )+b 


S (to) 
xx 


02.26) 


UT 


co‘ 


Once again, comparison of Equations (2.25b) and (2.26) shows 
that the non-linear device in cascade with an attenuation of 
0.5 will give a smaller noise output than the corresponding 
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linear element — - + , -U . regardless of the actual value of T , 

s m 

so long as Equation (2.24) is a reasonable approximation to 
Equation (2.23). 

Fig. 2.25 is a record of the square of the noise output 
(instantaneous noise power) of the linear and Clegg integrators 
produced by an input of zero-mean noise, and the results are 
in agreement with the above conclusions. 

Let us now consider the effect of zero-mean sensor-noise 
n(t) on the output c(t) (in the absence of command input), in 
the linear and non-linear systems shown in Figures 2.26 and 
2.27 respectively. 

We have , 


A S cc (w) 
N 1 (u) = 


S ( to) 
nn ' 


(b+1) 




linear system 
2 


P 1 (j«) 


'Xj 


x + (b+l)P 1 (jw) 
"jw 

^ |P 1 (jM) 1 2 


at high frequencies where 


P 1 ( jw) 


<< 1 , 


A S (u) 

M , % A cc 
N 5 (w) = j 

d S (w) i 

nn • 


j non-linear system 


(2.27) 
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(b+l)-e -J ’ wT 


m 


jw 


P 2 ( jw) 


1 + ( b + 1 - e ~ J ‘ a)Tm ) p 2 ^' (jj) 


? 2 

[b^+4 (b+ 1 ) sin (- 5 —)] ? 

% „ 5 |P 2 Uu)| 


to 


at high frequencies where 


P 2 ( jw) 




< < 1 . Thus , 


( 2 . 28 ) 


2 2 “^m 

N 0 ( to) [b +4(b+l) sin (— 5 —)] 

d % 

% 


N 1 (u) 


(b+ 1 ) ' 


P 2 (jco) 

P-j/ju)) 


(2.29) 


at high frequencies where 


p 1 (Jw) 




<< 1 , 


P 2 (jto) 

jw 


<< 1 . 


We thus find that 


N 0 (a>) I P 1 (j«) 

< 1 if | P P ( j u) ) 1 < 

N 1 (to) 1 2 


at high frequencies, 


i.e. if |L 2 (jw)| db < |L-,(jo))| dh - 6 , at high frequencies 


db 


where 


(b+l)P, (s) 

L ( s ) = 

x s 


the loop-transmission of the linear 
system of Fig. 2.26 


(b+l)P ? (s) 

L ( S ) = 

* s 


the linear loop-transmission corresponding 
to the non-linear system of Fig. 2.27. 


We shall show later that for given design specifications 
on the sensitivities of step- and disturbance-responses, the 
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use of the above non-linear device enables us to make a non- 
linear design in which the magnitude of (jtu) is at least 
6 db smaller, at high frequencies, than the magnitude of the 
loop-transmission L^joj) in the optimal linear design for the 
same problem. This will show that, for given specifications, 
the non-linear design achieves a smaller transmission of high 
frequency noise than is possible in the linear design. 



CHAPTER III 


PROCEDURE FOR NON-LINEAR DESIGN 


3 . 1 The Design Problem 

We now examine the application of the non-linear device 
of Fig. 2.14 in the design of a closed-loop system. We wish 
to develop a procedure for such design and to compare the 
results with the optimal linear design obtained by Horowitz 
and Sidi’s^^ method. 

The design problem is stated as follows : 

(a) The transfer-function P(s) of a linear, time- 
invariant minimum-phase plant is described by a set of para- 
meters a l5 ..., a^ (gain, poles and zeros). Each a^ (i = 1,... 
is known to lie within a given range < a i < B ^ (i = l,...p), 
but its precise value is not known. Thus P(s) is a particular, 
unknown member of a certain known family of functions. 

(b) Compensation must be designed in the two-degree-of- 
freedom structure of Fig. 3.1, such that, for all possible 
conditions of the plant P(s), the step-response of the system 
lies within certain specified bounds, such as, for example, the 
envelope shown in Fig. 3.2. 

(c) In addition, it is required that the response to a 
disturbance D at the output (Fig. 3.1) should satisfy certain 
damping requirements. 

(d) The magnitude of the loop-transmission in the linear 
design, and that of the corresponding linear loop-transmission 
in the non-linear design, should have an asymptotic slope of 
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( -20q) db /decade , as w tends, to infinity, where q is a given 
positive integer. In the linear design, q is the excess of 
poles over zeros in the loop-transmission. 

Also, the loop-transmission of the linear design and 
the loop-transmission corresponding to the non-linear design 
are both of Type-1. 

(e) The object is to arrive at a design which meets the 
specifications on step- and disturbance-responses, and, at 
the same time, minimizes the effect of feedback-sensor-noise, 
when there is no command input. A linear design which is 
quasi-optimal in this sense is furnished by the method of 
Horowitz and Sidi^^. 

It will be shown that, using the Clegg Integrator in 
parallel with a linear integrator, a non-linear design is 
possible which transmits less feedback noise to the output 
than the optimal linear design. 

3 . 2 Design Procedure 

The calculations that have been carried out for time- 
response with the non-linear device, and the results obtained, 
have been described in terms of the transfer- functions defined 
for the corresponding linear system. In keeping with this 
description of non-linear results in corresponding linear 
terms, we shall present the method of non-linear design as 
an adaptation of the method of linear design proposed by 
Horowitz and Sidi^ 1 ^. In fact, the method of non-linear 
design, for a given set of specifications, consists in obtaining 
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the corresponding linear transfer-functions as a linear design 
(by Horowitz and Sidi’s method) for a less restrictive set of 
specifications. The method of linear design is treated at 
length by Horowitz and Sidi in Reference [1], and, therefore, 
in order to describe the method of non-linear design, it will 
be enough to give an outline of their method and to indicate 
how it is adapted for non-linear design. 

3.2.1 Outline of Procedure for Linear Design^ 1 ^ 

The time-domain bounds, within which the step-response 

should remain for all possible plant conditions, are translated 

into corresponding bounds within which the function 

|T(jU)) | = | should remain for all plant conditions. 

1+L(jw) 

There is no rigorous method of obtaining such corresponding 
bounds, but Horowitz and Sidi propose a practical method of 
trial and error which produces bounds that lead to a successful 
design. 

The upper and lower bounds on |T(ju)|.^, in turn, lead to 

a bound A(w)db on the maximum variation that can occur, at 

fixed w, in | — 1 over all the possible plant conditions, 

1+L(jw) db 

since |T(jw)| db '= | F ( J co ) | ^ + • j kLsI“2 — | ^ and the pre-filter 

1+L(ju>) 

F(s) is fixed and independent of the plant. 

Therefore, at each w, the loop transfer-function L(Joj), 

considered at some chosen plant condition, should be such that 

the variation in I ^ . 4 — ) — | , due to changing plant 

1+L(ju>) 

conditions, should not exceed A(o>)db. The use of the Nichols 
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Chart is the obvious method of examining the variation in 


db 


due to variation in L(jto). Using templates of 

1+L(ju>) 

the variation of P(jto) over the plant conditions, and the 


bounds A ( <jo ) on the corresponding variation in 


L ( ;) to ) 
1+L (j to) 


db 


we obtain, for each to, a locus on the Nichols Chart as the 
lower bound for L(j<o). At high frequencies, however, the 


variation A(co) that is permissible in 


L(,joQ 


is greater 


1+L( jto) 

than the actual variation that occurs in L(jto), and, as a 
result, the step-response bounds impose no conditions on 
L(jto) at these frequencies, and there are no loci on the 
Nichols Chart. 

However, we must also take into account the specification 

that sets a limit to the overshoot in the disturbance response, 

determined by the transfer-function — 

2 1+L ( s ) 


. Using the 


(0 


second-order model 


n 


s 2 +26to n s+(o n ^ 


for 


L(s) 


, Horowitz 


1+L ( s ) 

and Sidi translate this into a lower bound on the damping 6 


in the second-order model, and, in turn, into a closed boundary 

on 1 — 1 , for all frequencies and plant conditions, based 

1+L ( j to) 

on the second-order relation 


Max 

to 


to 


n 


to 2 - to 2 + J26toto 
n. ° n 


26 


/l-6‘ 


The disturbance-response specification thus limits the maximum 


value that 


L(,j to) 


may assume, over all frequencies and plant 


1+L( jto) 

conditions; on the Nichols Chart this defines a- contour V (in 
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the notation of Horowitz and Sidi) whose effect is to limit 
the phase-lag that L(jw) can have until some frequency w at 

X 

which | L ( J oj )| is sufficiently small. The smaller the 

.X. 

permitted damping in the disturbance-response, the larger 
the phase-lag in L(jai) that is allowed by the contour V. 

Fig. 4.1 in Chapter IV shows the contours and V 2 corresponding, 
respectively, to 6 1 = 0.5 and 6 2 = 0*3, for a 40 db . variation 
in the gain of the plant. 

Thus, the specifications on step-response and disturbance- 

response, taken together, furnish a set of loci on the Nichols 

Chart as bounds on L(joj) at all w. A linear design that meets 

the specifications and is optimal, in the sense of minimizing 

k where k h - lim s q L(s), is one in which LCjco) lies on its 
s~*-°° 

corresponding locus at each w. The existence of such an 
optimal L(s) and the manner of its realization are shown in 
Horowitz and Sidi's paper 

3.2.2 Development of Procedure for Non-Linear Design 
Step-Response Specifications 

The step-response c(t) of the two-degree-of-freedom 

structure of Fig. 3.3, with the non-linear device (Clegg 

Integrator + — ) situated at the error junction of the feedback- 
. s 

loop is described by Equations (2.13) and (2.14) when we set 
R(s) = F 4- S i . It is clear that the non-linear response is 
identical with the corresponding linear response till t^, 
the instant of first reset of the. Clegg Integrator (in 
particular, if the Clegg Integrator is never reset, the linear 
and non-linear responses coincide for all t). This portion 
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of the non-linear response is thus determined by the 

corresponding linear transfer-function T(s) = and 

1+L( s ) 

therefore, its sensitivity to plant variations is, in fact, 
determined by the sensitivity of the linear transfer-function 
1 + l T s ) ( since p (s) is fixed and independent of the plant). 

We thus conclude that the variation of this portion, at least, 
of the non-linear step-response can be kept within the speci- 
fied bounds by designing the corresponding linear loop- 
transmission L(s) as in a purely linear design. 

As regards the response after the first reset, it is 
helpful to recall that a reset of the output y(t) of the Clegg 
Integrator at t i from y.^ = ytt^-), may be represented in the 
corresponding linear system by an impulse-input -y^SCt-t^). 

With the modified non-linear device (Clegg Integrator + 
situated at the error-j unction , the response, after t^ and 
until a possible second reset t 23 is found from Eqns . (2.13) 
and (2.14): 


c ( t ) = c T (t) - 

L 1+b 


h(t-t 1 ) , where 


r(t) = the output of the prefilter 
t 

c^(t) = | h(p)r(t-p)dp , the linear response. 


'o 


h(t) = impulse-response of 


L(s) 


1+L( s ) 


Yt = 


J 


e(t)dt = 
o J 


[r( t)-c( t) ]dt 


cn |o* 
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Since a reset occurs at t^, e(t) crosses zero at t-^ 
and changes sign. Consider the case when e(t) is positive 
for 0 < t < t^ and becomes negative after t-^, i.e. 

c(t) < r(t) for 0 < t < t^ and 
c(t) > r(t) for t-j, < t < t^ 

Then > 0. 

A[s n " q + ] 

Let L(s) = — 

[s n + ] 

where q is the excess of poles over zeros in L(s), and A > 0 
for minimum-phase L(s) . 

Then, from the initial value theorem of Laplace Transforms, 
we find 

h (l) (0) = lim s (l+1) = 0 , i = 0 ,...,( q-2 ) , 

s-* 00 1+L(s) 

and h (q_1) '(0) = A > 0. 

Hence, h(t) > 0 over an open interval (0,A) and thus, 

y-. 

c(t) = c,(t) - h(t-t,) < c L(t) 

1+b 

over an open interval (t-^, t-^+A) ; i.e. the reset at t^ tends 
to reduce the overshoot after t^. This is illustrated in Fig. 
3.4 . If a second reset occurs at t^, similar arguments show 
that it tends to reduce the undershoot that occurs after 
If there is no reset after t-^, then c(t) > r(t) for t > t-^, 

y-. 

c(t) - c T (t) = - h(t-t-,), t > t-. . Therefore, 

L 1+b l-i 



r 00 


[c(t)-c L (t) ]dt = - 


1+b 


h( t-t- L )dt 


*i_ 

l+b 


J o 


h(t)dt 


l+b 


lim 

s -*0 


L(s) 
1+L( s) 


l+b 


if L(s) Is of Type-1 or higher, showing that, on the average, 
the response c(t) has less overshoot than c^Ct) for t > t^. 

Thus, the effect of a reset is, in general, to reduce 
overshoot and undershoot in the non-linear response as 
compared to the corresponding linear response. In conse- 
quence, one expects the non-linear step-response. In 
general, to have smaller oscillations than the corresponding 
linear step-response, and, therefore, to remain within the 
same bounds as the corresponding linear step-response. 

An exception might occur for the plant condition in which 
the linear step-response corresponding to the non-linear 
design has near-critical damping. If the linear response 
has just enough overshoot to cause a reset in the non- 
linear response, the latter can have considerable undershoot , 
which might violate the lower bound of the step-response 
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specifications. Equation (2.19) gives an approximate* 
(second-order model) upper bound to such undershoot and 
enables b to be chosen so as to limit the undershoot to the 
desired value. 

A general conclusion from these considerations is that 
the non-linear response tends to be less oscillatory than 
the linear response for quite general inputs, exceptions 
occurring when a reset occasioned by an overshoot causes 
exessive undershoot or vice versa. Examples of linear and 
non-linear response to general inputs are shown in Figures 
4.11-a, b and c, for the linear and non-linear systems designed 
in Chapter IV. 


* Equation (2.19) applies to the step-response without a 
pre-filter, or, equivalently, to the disturbance-response; 
however, if the step-response of the pre-filter is much faster 
than that of the feedback-loop for near-critical damping, the 
expression is a good approximation even for step-response with 
the pre-filter. 

For instance, in the design example of Chapter IV, 


kL 2 

i+kLg has near-crit i cal damping for k = 1, and has then a 
bandwidth of about 4 rad/sec, while the pre-filter FpCs) = 


40 

s+4o 


has a bandwidth of 40 rad/sec. Figures 4.5-b and 4.6-b show 
that both the disturbance-response and step-response for k = 1 


have the same undershoot. 
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We thus conclude that sensitivity specifications on 
step-response can be met in the non-linear design by 
designing the corresponding linear loop-transmission for the 
same low-frequency boundaries on the Nichols Chart as in a 
linear design, with special consideration for the largest 
undershoot that can occur. 


Disturbance-Response Specifications 

The benefits of the Clegg Integrator over the linear 
integrator are secured, in the main, in the design to meet 
the disturbance response specifications. In linear design, 
the specification on the overshoot in disturbance response is 
translated into a minimum damping 5^ in a second-order model 
for the feedback-loop, and, in turn, into a boundary V-^ on 
the Nichols Chart, which limits, over a range of frequencies, 
the phase-lag allowable in the loop-transmission L-^(.jco) of 
the linear design. 

It is seen from Equations (2.17) and (2.18) that, for 
the same overshoot in the disturbance-response, the non-linear 
feedback-loop can have a smaller damping (in the corresponding 
linear transfer- function) than in the linear design. If 6 2 
is the damping corresponding to the non-linear design, equating 
the overshoots in the disturbance-responses of the linear and 


non-linear designs, we have 


7T 6 -j_ 

e ‘ 


b - 26 2 e 


V -1 

— " ^ 1 (tt— COS 5q) 

(1 -S§) 1/2 


1 + b 


TT<$ 2 

e ' 

(3.1) 
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which gives b as a function of 6^ and 5 2> 

Equation (2.19) shows that the choice of b also deter- 
mines the maximum undershoot that can occur in the disturbance- 
and step-responses. If this expression is used to set a lower 
bound to b, then Equation (3.1) may be used to find the 6 2 
that corresponds to a given 6^, and we note that S 2 < 6^. 

For a given 6^, the smaller the value of b, the smaller the 
value of <$ 2 . The non-linear design is then obtained by 
designing the corresponding linear loop-transmission L 2 (jw), 
using the boundary V 2 that corresponds to the 6 2 calculated 
from Equation (3.1). Since 6 2 < 6^, V 2 allows a faster 
reduction of |L 2 (jo>)| than is permissible in |L,(jw)|. 


Noise-Transmission 

It has been shown in Section 2.7 that if | L 2 ( J ai ) | < i|L 1 (ju) 

at high frequencies (i.e. at frequencies where 

— % L(juO), the non-linear design transmits less feed- 

1+L(jw) 

back sensor-noise than the linear design, if the significant 

part of the noise-spectrum occurs at these high frequencies. 

The use of the boundary V 2 in the non-linear design in place 

of the boundary V, in the linear design usually enables us to 

L 2 ( jw) ^ 

make the ratio lim [ 1 considerably smaller than — . 

or*- 00 L-^(jcj) 2 


3.2.3 Steps in Non-Linear Design 

The procedure for non-linear design may be summarized 


as follows: 



1 . 
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The specifications on step-response are translated into 
bounds on the variation, at different frequencies, of the 
same minimum-phase transfer-function T(s) that would be 
used in Horowitz and Sidi’s linear design. 

2. Using these bounds and templates of plant-variation at 
the different frequencies, boundaries are derived on the 
Nichols Chart for the loop-transmission as in the linear 
design. 

3. The damping 6^, which would be used to derive the high- 
frequency boundary in the linear design is calculated 
from the disturbance-response specification. Values for 
b and 6 2 are then obtained from Equations (2.19) and 
(3.1), the choice of b influenced by the largest under- 
shoot that can occur in the step and disturbance responses. 
The high-frequency boundary V 2 corresponding to 6 2 is 

then obtained. 

4. The linear loop-transmission L 2 (s) corresponding to the 

non-linear design is obtained as a linear design for the 
boundaries derived in Step-2 and the high-frequency 
boundary V 2 . , 

5. The series-compensation consists of the non-linear device 

(Clegg Integrator + |-) at the error junction, followed by 

s 

L (s) 

the linear transfer- function 2 , where 

P(s)[l + L ] 
s 

P(s) is the plant transfer function corresponding to the 
plant-condition used in deriving the boundaries in 
Step-2. 



The linear pre-filter F 2 (s) for the non-linear design 

is obtained from the equation 

F~( ju))L~( jw) 

T(jto) = 

1 + L 2 (j03) 

where T(s) is the transfer-function derived in Step-1. 



CHAPTER IV 


A DESIGN EXAMPLE 

4 . 1 The Problem 

The method of non-linear design is illustrated, and 
linear and non-linear designs are compared in the following 
example. 

Specifications : 

1c 

(i) The plant has the transfer-function P(s) = , 

s+4 

where 1 < k < 100. 

(ii) A pre-filter and series compensation are to be 

designed in order to keep the step-response of the two degree 

of- freedom system of Pig. 3.1 within certain bounds, for all 

values of k, 1 < k < 100. We assume here that such time- 

domain bounds lead to the following bounds on the variation 

of the magnitude of the minimum-phase trans fer- function 

T(s) = s that would be used In a linear design. 

1+L(s) 


A ( to) = Max |T(joj)| db - Min |T(ja))| db 


= Max: 


L(j«) 

1+LCjo) 


db 


Min 


L(ja)) 

1+L(ju) 


db 



48 . 
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TABLE 4.1 

Bounds On Variation in | T ( j oo ) | 


a) rad/sec 

Max|T(ju)) | db 

Min|T(j<o)| ab 

A ( to ) db 

1 • 

0 

-0.25 ■ 

o 

« 

[Y> 

ui 

2 

0.5 

-0.5 

1 

4 

1.0 

-4.0 

5 

6 

2.0 

-9.0 

11 

10 

2.0 

-18.0 

20 

20 

0.0 

-30.0 

30 

40 

0.0 

-40.0 

40 


(ill) The response to a constant output-disturbance 
should go to zero in the steady-state, with a maximum over- 
shoot of about 20$ of the constant disturbance. 

(iv) The loop-transmission L-^Cs) in the linear design, 
and the loop-transmission L 2 (s) corresponding to the non- 
linear design, should both be of Type-1 and should have an 
excess of four poles over zeros. 

4.2 Linear Design 

Since the plant variation is a variation of pure gain 
by a factor of 100, the corresponding template on the Nichols 
■Chart is a straight line of length 40 db , parallel to the 
magnitude-axis, for all frequencies. The bounds A(w) of 
Table 4.1 and the template of plant-variation together determine 
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the boundaries shown in Fig. 4.1 for L^(jw) at the different 
frequencies. k = 1 was used as the reference plant-condition 
in obtaining these boundaries. 

Corresponding to the 20 % overshoot allowed in the 
disturbance-response, 6^ % 0.5 and this leads to the high- 
frequency boundary shown in Fig. 4.1. Using these 
boundaries, the following loop-transmission L-^(s) is obtained 
for the linear design. 

2.3xl0 7 (s+ll) (s+35) 

L,(s) = 5 5 — (4.1) 

s(s+4)(s+30)(s+180)(s+216s+360j 


Fig. 4.1 shows that L^(jo)) lies quite close to its 

boundaries, and therefore, it is close to the optimal design 

in the sense of Horowitz and Sidi. 
kL,(jo>) 

With known, for 1 < k < 100, the following 

l+kL^jw) 

pre-filter F^(s) for the linear design is obtained with the 
help of Table 4.1 on pg. 49. 


F 1 (s) 


70 

s + 70 


(4.2) 


This completes the linear design, and the Bode plot of 
L^Cs) is given in Fig. 4.2. 

4 . 3 Non-Linear Design 

The loop-transmission LjCjw) that corresponds to the 
non-linear design uses the same low-frequency boundaries as 
were obtained for L-^(jw) (with the same reference 
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plant-condition). When <$ 2 = 0*3 is tried. Equation (3.1) 
gives b = 2.7 (the calculation is shown in the Appendix), 
which, from Equation (2.19), would limit the maximum under- 
shoot to about 20$. The high-frequency boundary V 2 
corresponding to 6 ^ - 0.3 is shown in Fig. 4.1. 

L 2 (s) is now obtained as a linear design with V 2 
replacing as the high-frequency boundary, and is given by 

4.6x10 6 (s+18)(s+45) , h 

L (s) = — 

s(s+4) (s+50) (s+8o) ( s^+208s+260^) 


kLp ( j w) 

With known, for 1 < k < 100, the pre-filter 

l+kL 2 (ju>) 

F 2 (s) for the non-linear design is chosen so as to make 
F 2 ( jw)kL 2 (jm) 

vary, as k varies from 1 to 100, between the 

l+kL 2 ( Joj) 

limits given in Table 4.1 on pg. 49, over the range of fre- 
quencies listed in the Table. We choose 


F,(s) = 


s+4o 


2 7 

The non-linear device (Clegg Integrator + — —) is 

s 

situated at the error-j unction of the feedback loop, followed 
by G(s), the linear part of the compensation and the plant 
P(s) . G(s) is determined by the equation 


L 9 (s) = P(s)G(s) t.3± b l = P(s)G(s) ^1 
^ s s 


This completes the non-linear design and the Bode plot of 
L 5 (s) is given in Fig. 4.2. The Bode plots of T-, = F,(.s) 


1+kL (s) 
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kLp(s) 

and T ? (s) = F p (s) , for k = 1 and k = 100, appear 

_ l+kLp(s) 

in Figs. 4.3 and 4.4, respectively. 

4 . 4 Results 

Figures 4.5 and 4.6 show the step- and disturbance-responses 
of the linear and non-linear designs, for several values of k, 
as well as the responses obtained when Lp(s) is used as a 
linear design. 

We note that the Lp-non-linear responses fall within the 

bounds of variation of the L^-linear responses and meet the 

same specifications as the linear responses; the Lp-linear 

responses , however , violate those bounds. 

The transmission of feedback sensor-noise to the output 

in the linear and non-linear designs is compared in Fig. 4.7, 

which is a record of the square of the actual noise output 

produced by zero-mean sensor noise (with a bandwidth of about 

lOOLp ( s ) 

fifteen times that of ) , and is thus a record of 

1+100L 2 (s) 

instantaneous noise-power at the output. The noise at the 
output is seen to be much smaller in the non-linear design 
than in the linear design, and as a result, saturation of 
forward elements due to sensor-noise is much less of a problem 
in the non-linear design than in the optimal linear design. 

Figures 4.8, 4.9 and 4.10 show linear and non-linear 
responses to sinusoidal and ramp inputs and to a combination 
of command- and disturbance-inputs. Figures 4.11-a,b,c compare 
linear and non-linear responses to some other input waveforms. 

Though a Type-0 plant and pure gain-variation were 
considered in this example, the derivation of the procedure 
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for non-linear design has been quite general, with just one 
assumption made — that the device (Clegg Integrator + ^-) is 
situated at the error-j unction ; however, this imposes no 
restrictions on the design, and we see that it is also 
applicable to other kinds of plants and parameter-variations. 

4 • 5 Conclusion 

It has been shown that the non-linear device consisting 
of the Clegg Integrator in parallel with the linear integrator 
is a useful element to include in the compensation for a 
linear, minimum-phase plant with uncertain parameters, to 
achieve specified bounds on its response. For given plant- 
ignorance and bounds on response, it has been shown that, 
using this non-linear device, compensation can be designed in 
a two-degree-of-freedom structure in which there is less 
transmission of feedback sensor-noise to the plant and its 
output than is possible to achieve in a linear design. A 
method has been obtained which enables the non-linear design 
to be carried out in terms of an equivalent linear design by 
the method of Horowitz and Sidi. 
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r(t) 



L(s) = P(s)G(s) 


T(s) = 

1+L(s) 


B ( s ) = - — 


J_ 

PCs) 


L(s) 
1+L ( s ) 


Fig. 1,1: A Two-Degree of Freedom Feedback System 






:aj^gSSaBBBM^ ggg^ a^fegji«a- 5tewga a ggP M^ ^5S^5SSB^§ a£g-^Sg BgB &^.i 

gg g s u ai«^ 1 1 L P t i ■ « ■ BMB mffiaaMjBHMgB gg »^hwi«ii^i i ii i h iii i h |gg ttgBSB 31 

al|SnHSilsssi 9 H 

!SaffiaaEKga!g™iSigSBB3K™~”g'~»BB!^JS:gB^tj«ggag aa8M p F ^- g; ‘ :;g> * <M * g * gg;< ^ 1 **- 1 " 

sssSSl 


2 =I=iIIi=iI=i 


gi 




:! i:lf i: f h 



■■■■ 

Ill 











::S±t 


ssais 










> ■■■■■■■■■■■■■■■ 
!■■■■■■■■■■■■■■ 

.■■■■■■■■■■■■■■■■■ 

■Hummummimi 



SKI 

sss 



=asa i 

a£S 


Eli 














raw 


s»: 






Il!liiai|H|liII 

■ssssiKKKKKKIi 





!iiH 


¥f 


gp -T 

aa 33 ^;i 












■fcssss^sssfl 

EssKSiEfficSisSi 

5 555555 j a *.a^ 3 MI 

:==:== sHaslil 







m 






— i 

h t — i 

r- : ” 

! 

3 j 

it. 


9 

























57 



o 

C\J 


o 


o 











Compute 



Simulation of the Clegg Integrator 
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R(s) A(s) R(s) A(s) 


Fig. 2.10: Linear Equivalent of System of Fig. 2.9 





















































Bandwidth of noise % 10 rad/sec 
Peak-to-Peak value % 160 volts 




Fig. 2.25: Output Noise-Power of Linear and Clegg 

Integrators with Zero-Mean Input Noise 
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25,50,75,100: 1 cm = 0.05 sec. 
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Pig. 4.7: Output Noise-Power of i+iq'OL 
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Lnusoldal Input , f = 20 c/s 

rdinate Scale: 2 volts/line 

ime Scale: 1 cm = 0.1 sec. 

esDonse Scales (Same for all three 
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Fig. ^.8: Sinusoidal Resp 
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Scales (same for all four figures) 
Ordinate: 2 volts/line 

Time: 1 cm = 0.1 sec. 



T 2 - Non-Linear 



T 2 - Linear 


Fie. 4.11-a: Linear and Non-Linear Responses 
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Ordinate: 2 volts/line 



Time: 1 cm = 0.1 sec. 

Response Scales (for all three 

responses ) 

Ordinate: 0.5 volt/line 

Time: 1- cm = 0.1 sec. 


T^ - Non-Linear 



Fig. *J.ll-c: Linear and Non-Linear Responses 
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APPENDIX 


3M n 

A.l Calculation of — — - (6 b) 

3 6 


From Equation (2.17) on page 24, setting <j> = cos ^6 , 
M n (6 ,b) = M n (cos<t> J b) 

b tt cot <|) 2cos(j) e -<j>cot<f> 



1+b 

e 

- e 

1+b 

9M 

n 

b 

(-rrcosec 

2 A -7rcotcf>s 2e~^ cot<t> t , xtJ . x 2, 

<J>e T ) - {-sin<j>+<f>cos(f>cosec <j> 

3 <f> 

1+b 


1+b 




-COS<J)COt(j>} 



1 

o 

rTrbe -7TC0 ^^+2e - ^ G0t ^ ( sine})- cf> cos <f>)l 


(1+b ) sin <j> 

L J 


> o 

, o < 

4> < - (i.e. 0 < 6 < 1) . 

2 - 

Now , 




3M n 

_ 3M n 

3 4> 


36 

3 tj) 

36 



; 3M 

n 

3 ( cos -1 

6) 


3<f> 

36 



1 

3M n 



- 

< 0 , 0 < 6 < 1 . 


/l-6 

2 3(f) 



A. 2 Calculation of b in Non-Linear Design of Chapter IV 

We use Equation (3.1) on page 44 with <J>^ = cos -1 6 1} 
. A -1. 

<p g = cos 6 £ • 
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r i 


e “TT cot <{> x = _1_ U _ 25 e cot<*> 2 (ir-<t> 2 ) 

1+b l: 2 


; -ircot<|>2 


1 +b 


be -frcot^ 2 _ 26 ^ e -<p 2 Cot(P 2 


Prom the linear design, S JL = 0.5 and therefore, «(, = 60° and 

-ircotij)-! ■ c , 

0.164. With & 2 = 0.3 In the non-linear design, 

<j>2 = 72.5° = 1.265 radians, and e “ Trcot< ^2 = 0.374. 
2S 2 e-*2 c ° t *Z = 0.403.' Thus , 

0.164 = — [0.374b - 0.403] 

1+b 

b = 0^03 + 0.164 m 0,567 = 2 

0.374 - 0.164 0.21 



